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We investigate the entanglement-related features of the eigenstates of two exactly soluble atomic 
models: a one-dimensional three-electron Moshinsky model, and a three-dimensional two-electron 
Moshinsky system in an external uniform magnetic field. We analytically compute the amount of 
entanglement exhibited by the wavefunctions corresponding to the ground, first and second excited 
states of the three-electron model. We found that the amount of entanglement of the system tends 
to increase with energy, and in the case of excited states we found a finite amount of entanglement 
in the limit of vanishing interaction. We also analyze the entanglement properties of the ground and 
first few excited states of the two-electron Moshinsky model in the presence of a magnetic field. The 
dependence of the eigenstates' entanglement on the energy, as well as its behaviour in the regime of 
vanishing interaction, are similar to those observed in the three-electron system. On the other hand, 
th e entanglement exhibits a monotonically decreasing behavior with the strength of the external 
magnetic field. For strong magnetic fields the entanglement approaches a finite asymptotic value 
that depends on the interaction strength. For both systems studied here we consider a perturbative 
approach in order to shed some light on the entanglement's dependence on energy and also to clarify 
the finite entanglement exhibited by excited states in the limit of weak interactions. As far as we 
know, this is the first work that provides analytical and exact results for the entanglement properties 
of a three-electron model. 



I. INTRODUCTION 

Entanglement is an essential ingredient of the quantum mechanical description of Nature [H-O- Besides its cen- 
tral role for the basic understanding of the quantum world, entanglement constitutes a physical resource admitting 
numerous technological applications. The study of entanglement sheds new light on the mechanisms behind the 
quantuni-to-classical transition [4| as well as on the foundations of statistical mechanics Q . On the other hand, the 
controlled manipulation of entangled states of multipartite systems is fundamental for the implementation of quantum 
information processes, such as quantum computation @, Quantum entanglement is also relevant in connection 
with the physical characterization of atoms and molecules. The exploration of the entanglement features exhibited 
by atoms and molecules is a captivating field of enquiry because these composite quantum objects play a central role 
in our understanding of both Nature and technology. In point of fact, the entanglement properties of atomic systems 
have been the subject of considerable research activity in recent years 0, [sl-flTj. This line of research is contained 
within the more general one aimed at the application of information-theoretic concepts and methods to the study of 
atomic and molecular systems [isl - lsoj . 

Some of the most detailed results on the entanglement properties of atomic systems, particularly in the case of 
excited states, have been obtained from analytical investigations of soluble two-electron models @, [l^. Partial results 
were also obtained numerically for the eigenstates of helium-like systems, employing high quality wave functions 
[isj . Some general trends are beginning to emerge from these investigations. It is observed that the amount of 
entanglement of the atomic eigenstates tends to increase with the concomitant energy. It also increases with the 
strength of the interaction between the constituent particles. On the other hand, the entanglement of excited states 
shows an apparent discontinuous behaviour: it does not necessarily vanish in the limit of very small interactions 
It would be desirable to extend these studies to more general scenarios, particularly to models consisting of more than 
two electrons, or involving magnetic fields. The aim of the present contribution is to investigate the entanglement 
properties of the eigenstates of the exactly soluble Moshinsky model [3l| , extending previous works to the cases of a 
three-electron system and a three-dimensional two-electron system in a uniform external magnetic field. 

The paper is organized as follows. In section 2 we briefly discuss entanglement in systems of identical fermions. 
We review the measure used in order of quantify the amount of entanglement of pure states, focusing on appropriate 
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measures for two- and three-electrons systems. In section 3 we investigate the entanglement properties of the eigen- 
states of the Moshinsky model with three electrons. The entanglement features of the three-dimensional Moshinsky 
model with two electrons in the presence of a uniform magnetic field are studied in section 4. Then, in section 5 we 
consider a perturbative approach to clarify some entanglement features found in the previous models. Finally, some 
conclusions are drawn in section 6. 



II. ENTANGLEMENT MEASURE 



Correlations between two identical fermions that are only due to the antisymmetric nature of the two-particle state 
do not contribute to the state's entanglement [32- 37]. The entanglement of the two-fermion state is given by the 
quantum correlations existing on top of these minimum ones. A practical quantitative measure for the amount of 
entanglement exhibited by a pure state ji/)) of a system of N identical fermions is (see [sl] and references therein) that 
given (up to an appropriate multiplicative and additive constant) by the linear entropy of the single particle reduced 
density matrix p^, 

e{m^l-NTr[p% (1) 

Notice that, according to this entanglement measure, a pure state that takes the form of a single Slater determinant 
has no entanglement. The measure ^ is normalized to adopt values in the interval [0, 1]. 

We shall apply the measure given by Eq. ([T]) to a pure state |^') of a one dimensional system consisting of three 
spin-i fermions (electrons). This pure state has an associated wave function given, in self-explanatory notation, 
by *(a:;icri,a;2(T2, xaCTa) = {xiai,X2cr2,X3a3\^) with \xiai,X2a2,X3a3) = jxi, 2:2, xa) (g) |cri, 0-2, 0-3). Here 2:1,2,3 are 
the coordinates of the three electrons and the dichotomic variables cri.2,3 (each adopting the possible values ± and 
corresponding to the Sz component of spin) describe the spin degrees of freedom of the three electrons. In order to 
evaluate the amount of entanglement of the system we have to compute the following integrals, 



{Xiai\pr\x[a[) = ^ / {xiai,X2(T2,X3a3\p\x[cr[,X2<J2,X3<T3)dx2dX3 (2) 

where {xi<7i\pr\x[a[) are the elements of the one-particle reduced density matrix pr{xi, x'l), p = and 

{xi(Ji,X2(J2,X3,az\p\x'i(7[,X2(72,Xza3) = ^!{xiai,X2(72,Xz(ys.)^*[x'^(7'i,X2(J2,X:>,(7z)- 

The square reduced spin density matrix is given by 

/oo 
(xifTi \pr\xa) {xa\pr \x[a[)dx 

and finally the expression for the trace is 

/oo 
{xa\pl\xa)dx 

o=nz -°° 

In the three-electron case it is not possible to find totally antisymmetric factorizable between coordinates and spin 
wave functions, however in the two-electron case, following [8] we focus on states with factorized wave functions. The 
corresponding density matrix takes the form 



p = p(-) ® p(^) (3) 
and then, the entanglement measure evaluated on these states is given by 

e = l-2Tr[{pi^yf]Tr[{p(/yf], (4) 

where pi'^^ and pi'^^ are the single-particle reduced coordinate and spin density matrices. So, in the case of two-electron 
system studied in section 4, we consider separately the cases of parallel and antiparallel spin wave function. In the 
case of parallel spins, described by | -I- or | ), the coordinate wave function must be antisymmetric and we have 
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Tr[(pr )^] = 1- On the other hand if we have antiparallel spins, we can distinguish two cases: symmetric coordinate 
wave function with spin wave function of the form -^(1 H — ) ~ I — I")) oi antisymmetric coordinate wave function 

with spin wave function -^(| H — ) + | — h)), both of them with rr[(pf'*-')^] = i. And finally, to calculate the amount 
of entanglement we will compute the integrals 

(ri|p(='|r'i) = / (rir2|p(^)|r'ir2)dr2 = / vl/(ri, ra)** (r'l, r2)dr2 (5) 

and the trace of the coordinate part is 




III. THE THREE-ELECTRON MOSHINSKY ATOM 

The Moshinsky atom [sT] is a system formed by harmonically interacting particles confined in a common, external 
isotropic harmonic potential. The total Hamiltonian of the one-dimensional Moshinsky atom with three electrons is 

1 / d"^ d"^ \ 1 1 

^ " " 2 [d4 ^ 9^ ^ dxl) ^ 2^'^^' + •''2 + ^ 2^''li^i - + {X2 - x:if + (xg - x^f] (7) 

where xi, X2 and are the coordinates of the three particles, uj is the natural frequency of the external harmonic 
field, and A is the natural frequency of the interaction harmonic field. The positive sign in the last term describes an 
attractive interaction between the electrons and the negative a repulsive interaction. We use atomic units (me = 1, 
h= 1) throughout the paper, unless indicated otherwise. 
Introducing the Jacobi coordinates for three particles, 

Ri = ^{xi+X2+X:i), R2^ ^{-2xi+X2 + Xz) and i?3 = -^(x2 - xa) (8) 

the Hamiltonian separates in the following way, 

where jSi — uP' and /32 = /^s = = ± 3A^ (again, the + sign corresponds to an attractive interaction, while the 
— sign corresponds to a repulsive one). In the case of a repulsive interaction it is necessary to impose the constraint 
A < in order to obtain bound eigenstates. The general eigenfunctions of the system are 

«'(X1, X2, Xa) = 4'(i?l, i?2, Rz) = (^l)*n„, (i?2)*„H3 (^s) (^0) 

with 

(^0 = ^ , (11) 

where Unix) denote the Hermite polynomials. The eigenenergies of these states are 

E = Er-^ + Er,^ + Er., = UJ (^R^ + ^) + ^ ("-^2 + ^^^.3 + 1) (12) 

We will denote by \nR-^nR^nR^) the eigenstates of the Hamiltonian ([9|), which are characterized by the three 
quantum numbers n/j^, ur^ and ur^. To fully define the three-electron system's eigenstates we must take into account 
combinations of such functions of the coordinates together with the spin ones |(TiCT2C3) to obtain total antisymmetric 
wave functions. In this case the wave functions corresponding to the energy eigenstates cannot always be chosen 
to be separable between coordinates and spin and there are no spin functions totally antisymmetric by themselves. 
The Hamiltonian commutes with the spin observables, since it does not explicitly involve the spins. In particular, it 
commutes with the total z-component of spin angular momentum Sz- Consequently, it is possible to choose energy 
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eigenstates that are also eigenstates of Sz- It is plain that the wave functions associated with these eigenstates can 
always be written (up to a global normalization constant) in one of the forms 

|$(++-'>| + +-) + |$(+-+)>| + -+> + |$(-++)>| - ++), (13) 

|$(+++))| + ++), (14) 

or in the forms obtained substituting + by — (and viceversa) in the above expressions. In (|13m4[) the 
kets |<I)('^i''^2,o-3)^ correspond to the translational degrees of freedom and have associated coordinate wave func- 
tions ^'•'^^•''''•''■'\xi,X2,X3) = (a;i,a;2,a;3|$('^i''^='^'^2'). For the states (jlSIIUp to be fully antisymmetric the co- 
ordinate wave functions $'^°'i^°'^'°'-'^(xi,a;2,a;3) must satisfy the following set of relations. If ai = (72 we 
must have $^'^i''^^''^^''(a;2, xi, 2:3) — —^^'^^'"'^'''^\xi,X2tX3), (that is, in this case the coordinate wave function 
<j)(cri,CT2,CT3)('2.^^ 2:2, X3) has to be antisymmetric with respect to xi and 2:2). On the other hand, if cti — —02 we 
must have, $^'^i''^^''^^^(a;2, a;i, 0:3) = — <I>^°'^''^i''^^^(xi, a;2, xs). Similar relations must hold in connection with the pairs 
of labels ((72,(73) and (f73,cri). These relations imply, in particular, that the wave function $("'""'""'"''(a;i,X2,a;3) (and 
also $*■ ■*(a;i,a:2,X3)) must be fully antisymmetric in the three coordinates X\,X2tX3. Finally, it is clear that in 
order to be energy eigenstates the states p3m4p must involve spatial wave functions (xi, X2, 0:3) that are 

themselves eigenfunctions of the Hamiltonian ([7]). In particular, the three coordinate eigenfunctions associated with 
(jlSp must be eigenfunctions of ([7]) corresponding to the same energy eigenvalue. The ground state and few excited 
eigenstates of the three-electron system that we are going to study in the present work do not correspond to the form 
P^ . Thus, we are going to restrict our considerations to eigenstates of the form A direct way to construct the 
ground and first few excited states according to the structure (|13p is to use combinations of the forms, 



+ {\nR'^nji>^nii0 - \nji>^nji,m^ji>^)^ \ - ++) (15) 

or 

\nin2n3)=jV' \nR'^nii,^nji,^)\ + +-) + \nji,^nji,^nR0\ -\ + InR^nR^nR^)] - ++) , (16) 

where 

Ri = R[ = R'i, R'2 = -^(xi - 2x2 + X3), R'i = -^(xi + X2 - 2x3), 
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= — (x3-xi),i?^' = — (X1-X2), (17) 



ni = riR^ = ur:^ = nR:> , ^2 = ur^ ur/^ = ur'^ , 713 = ur^ = ur^^ = ur'^ , 

and A/", M' are appropriate normalization constants. Note that the three spatial wave functions corresponding respec- 
tively to the three kets \nR^nR^nR^) , |?^_R;n./^^ and \nRnnRqnRn) (which appear in (fT5|) and in ([T6|l ) are obtained 
via cyclic permutations of the particles coordinates xi, X2, X3 in the definition of the Jacobi coordinates. Therefore, it 
is evident that these three spatial wave functions are eigenfunctions of ([7]) sharing the same eigenenergy. We will use 
combinations of type (jlSp if the quantum number 71,3 is even, and of type (|16p when it is odd, ensuring in this way the 
antisymmetry of the wave function. As already mentioned, we have chosen these states because they are also special 
in the sense that they all are eigenstates of Sz ■ States of the forms and correspond to a wave function with 
total spin Sz — si^^ + si'^^ + 51'''' ~ +^ oi the three-electron system but one can also construct eigenstates of the 
same type with total spin Sz ~ As the entanglement of the Sz — —-^ states is the same as the entanglement of 
states with 6*2 = -I- i, in the rest of this work we will mainly focus on states with Sz = +-^- 

We must remember that these states are written in Jacobi relative coordinates of a three-particle system and 
the quantum numbers rii, 712 and refer to these coordinates. However, to determine the amount of entanglement 
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between the particles we have to express the wave functions associated with the eigenstates in terms of the coordinates 
and spins of the particles, 

'^mn27i3{^i'^i,X2cr2,X3a3) = (a;iai,a:2cr2, 130-3 |nin2n3) (18) 

In the case of the eigenfunctions (US]) of the Moshinsky system the entanglement measure e can be computed in 
an exact analytical way. However, for highly excited states the corresponding expressions become very awkward. 
Therefore, we are going to calculate this quantity only for the ground state and the first and second excited states. In 
each case we shall provide the final closed expressions for e (arising from the evaluation of the aforementioned integrals) 
and discuss the behavior of the eigenstate's entanglement. The value of s corresponding to the state \n1n2n3) (with 
Sz — +^) will be denoted by en^njng. As a compact alternative notation for the alluded state we shall also use 

hl"2»^3)flii?,2-R3- 

We compute the state's entanglement in terms of the dimensionless parameter t = ^, which constitutes a measure 
of the relative strength of the interaction between two particles in the Moshinsky system. Remark that the system is 
decoupled when r = 0. The larger the value of r, the larger is the (relative) contribution of the interaction term in 
the Moshinsky atom. 



A. Ground state |010)rj^r2R3 



Let A — vl±3r^, and using the right spin combination given by (|15p we can express the entanglement of the 
ground state in terms of the parameter r as 



£010 — 1 



y/2A + 5^2 + 2^43 
4(2 + 5^ + 2^2)3 



59 + 232A + 390^2 + 232^^ + 59^" 



(19) 



We see from (|19p that the entanglement of the ground state depends upon the parameters of the Moshinsky atom 
only through the dimensionless quantity r. Decoupling the system, that is, making r — )■ (which corresponds, for 
instance, to A — or equivalently A ^ ui) makes epio = showing that in the decoupled system the ground state is 
not entangled. On the other hand, with maximum coupling t — ^ 00 (t — > for attractive (repulsive) interactions 
we find that £010 = Ij that is, the entanglement measure adopts its maximum possible value. 



B. First excited states |110)r^r2R3 and |011)r^r2R3 

The first excited state in energy, when the system is coupled (r > 0) and with attractive interaction, is \iW) n-^n^Ri 
and the next one with higher energy is R-^n^fi^ , the excitation order is reversed in the case of repulsive interaction 
Both states have the same energy when we decouple the system, that is, when r — >■ 0. For these states, using 
and (|16p respectively, we have 



£110 — 1 



^1/2 



4 (2 + A)^/^ (1 + 2A) 



9/2 ■ 



'177 + 1034A + 6213^2 + 12582^=^ + 15392^" + 12582^5 + 6213^^ + 1034A^ + 177 A^ 



and 



^011 — 1 



640 (2 + Af^^ (1 + 2A) 



x9/2 



(20) 



(3057 + 24608^ + 93180^^ + 196704^^ + 251366^" + 196704^-^ + 93180^^ + 24608A^ + 3057^4^ 



(21) 



Decoupling the system makes £011 = £110 ~ 27' limit of a decoupled system the first excited states are 

entangled. On the other hand, with maximum coupling we find that eon ~ £110 — 1, that is, the entanglement is 
maximum. 
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To both states {\011)riR2R3 and \110) r^r^Rs) having Sz = +5, which we wih denote by |011)+ and |110) + , one can 
associate the states |011)_ and |110)_ respectively with the same energy and same entanglement but with Sz = — ^• 
Then, as these are degenerate states because the energy does not depend on the spin, we compute the amount of 
entanglement of a combination of them in the following way: 

1*011) = cos6l|011)+ + sine'lOll)- 



1*110) = cos0|llO)+ + sin6'|110)_ (22) 

where < 9 < 2tt. These states exhibit an amount of entanglement that is independent of the parameter 6. To 
understand this behaviour let us consider the unitary transformation U (acting on the single-particle Hilbert space) 
defined by, 

U\^k)\+) = \<Pk)\p), fc = l,2,... 

Um\-) = \c^k)\n), fc = l,2,..., (23) 

where (|efc) = |(/)fc)|±), fc = 1, 2, . . .) is a single-particle orthonormal basis (with the kets \(pk) corresponding to the 
spatial degrees of freedom) and 

\p) = cos 6*1 -f) - sin 6*1-), 

|n) = sin6'|-F) -hcos6'|-). (24) 

It can be verified after some algebra that, 

1*011) - (u<s>U(E,u)\on)+, 

1*110) = (U(E)U(E)U^\110) + . (25) 

Now, it is clear that the amount of entanglement of a three-fermions state does not change under the effect of unitary 
transformations of the form U ® U ® U and, consequently, the entanglement of the states defined in ((22)) does not 
depend upon 9. 



C. Second excited states |210)riR2R3, |111)r^r2R3, |012)r^r2R3, |021)riR2R3 and |003>RiR2R3 

For these states we have that the lowest-energy second excited state when the system is coupled and with attractive 
interaction is \2\Q) r^r^r^, the next one with higher energy is \\\\) r^r^r^, and the following three states, all of them 
with the same energy, are \Q12) r^r^r^, \{)21) r^r^r^ and r^r^^r^. All these states have the same energy when the 
system is decoupled. 

Defining the parameter B = (^2+A)^^r^(i+2A)^'^'^ ' "^"^S C3 states 1210)^1^2^3 and \012) r^r^^r^, and 

(flG]) for the states r^r^r.^, \Q21)r^r2r^ and r^r^r^, we found that 

£210 = 1 (2419 + 19480^ -f 218138/12 + 564200^1^ 4- 1466241^"* + 2943840^'^ -I- 3743124^^-h 

16 

-^-2943840^^ + 1466241A® + 564200yl^ + 218138A^" -I- 19480A" + 2419^^^) , 

£111 = 1 - — (9171 + 80546A + TOOSSSA^ + 2659770^3 -t- 6668841A* + 11416740^^ + 13615794^^-^ 
64 ^ 

-1-11416740^^ -f 6668841A® -I- 2659770^^ -I- 700555^^° + 80546A" + 9171^^^) , 

£012 = 1 (42739 + 506008^ 4- 3123242^2 + 11179160^^ + 26922957^^ + 44982480^^ + 53234988/1^-^ 

256 

+44982480A^ -I- 26922957A*^ + 11179160/1^ + 3123242^^" -I- 506008/1" + 42739/1^^) , 
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£021 — 1 — 



B 
4096 



(727363 + 8982520A + 54219206^^ + 196856600^3 + 469858317^1* + 776694000^^+ 



-915625428^*^ + 776694000^^ + 469858317^* + 196856600^^ + 54219206^^° + 8982520A" + 727363^^^) 



and 



^003 — 1 



B 
4096 



(762395 + 9419160^ + 61156086^^ + 232139320^^ + 576896949^'* + 982782000^^ 



+1171448436A^ + 982782000^^ + 576896949^^ + 232139320^^ + 61156086^^° + 9419160^" + 762395^^^) , 

4 43 1 

Taking the limit for the decoupHng case of the system, makes em = £210 = £012 = 7:, £021 = 77^ and eoo3 = 7 

9 ^ 108 ^ 4 

showing again that these excited states are entangled in the decoupled system. In the maximum couphng hmit we 
find for all second-excited states that the entanglement reaches again its maximum value. 



c 
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FIG. 1: (Color online) Entanglement of the ground, first and second excited states of one-dimensional Moshinsky atom with 
three electrons attractively interacting. All depicted quantities are dimensionless. 



The behaviour of the eigenstates' entanglement as a function of the parameter r (which corresponds to the relative 
strength of the interaction between the two particles) is depicted in Figure 1 for an attractive interaction and in 
Figure 2 for a repulsive interaction. 

Comparing Figures 1 and 2 one observes that in the repulsive case (Fig. 2) maximum entanglement is reached when 
the parameter r approaches the finite limit value ^ 0.577. In the attractive case (Fig.l) entanglement behaves in 
a different way: maximum entanglement corresponds to the limit r — > 00. This difference between the attractive and 
the repulsive cases is due to the fact that the Moshinsky model with repulsive interaction admits bound states only 
for T-values in the finite range [0, ^). On the other hand, in the attractive case the Moshinsky model admits bound 
states for all r > 0. In the case of the repulsive interaction the eigenstates of the system are no longer bounded for 
T > Thus, the eigenstates exhibit a qualitative structural change at the "critical" value Tc = resembling 
a quantum phase transition. A similar situation occurs in the case of the Moshinsky atom with two electrons in a 
uniform magnetic field studied in the next section (see Figure 4). This system, when the interaction is repulsive, 
admits bound states only for r- values smaller than the critical value Tc = 1. 
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FIG. 2: (Color online) Entanglement of the ground, first and second excited states of one-dimensional Moshinsky atom with 
three interacting electrons for the repulsive case. All depicted quantities are dimensionless. 



IV. TWO-ELECTRON MOSHINSKY ATOM IN A MAGNETIC FIELD 



The Hamiltonian of the two electron Moshinsky atom in a three dimensional space is 

where subscripts 1 and 2 denote each of the electrons. As before, the positive (negative) sign refers to a attractive 
(repulsive) interaction between the electrons. To study the presence of a uniform magnetic field acting on the system, 
we perform the following change in the Hamiltonian 

6 6 1 

Pi ^ Pi + - A and P2 P2 + - A with A = - (B A r) 
c c 2 

being B the magnetic field. Assuming that the magnetic field is homogeneous and have z-axis direction, that is 
B = Bz, we can write: 

2 

eS \ , n eB , 



Pi -^PX+ y^cj +yi) + — {xiPyi - ViPxi) with z = 1, 2 (26) 
By replacing (j26p in the Hamiltonian Hm and setting atomic units {m^ = h = 1, c — 1/ a), we obtain 

H = \{pI +pI) + + rl) +^^{xl+yl + xl+ yl) + 6(Li, + L2.) ± ^(^i - ^2? (27) 

where 

6 = — ; = [xiPyi - yiPxi) and = {xi,yi,Zi) with i = 1,2 (28) 
We change the variables to the center of mass (CM) and relative coordinates, i.e. 

R=i=(ri+r2) y r = -^(n - r2) (29) 
respectively. This transformation satisfies the relations 

p\+P2=Pr+pI and Li^ + L2Z = Lr^ + Lrz ^ {RxPB.y - RyPB.x) + {rxPry - TyPrx), (30) 



and therefore, introducing equations (l29l) and ([30]) in the Hamiltonian ((27|) we obtain 

H = \{pI+pI) + + r') + y (i?^ + Rl + rl+ r^) + 6(Lfl. + L„) ± yr^, 

which is separable in the CM and relative coordinates, so that we can express as 

H = Hfj + Hj, 

where 

Hr = ^(pL + ply) + ^^^^{Rl + Rl) + \{pI. + ^^Rl) + bin, 

and 

^^r- = +pIv) + ^^^^(^' + + \^pI^- + ^'^') + ^^rz ± 

Introducing the following dilation canonical transformation for the Hamiltonian H^i 

p'ji,^{uj^ + b'^)-ipR,,Ri^{uj'^ + b'^)iR, with i^x,y 



and 



we obtain 



p'rz=^ ^PRz, K=uj^R^, 
p^ = (w^ + 6^ ± A^)~3pH, = (w^ + 6^ ± A^)3ri with i = x,y 



CO Z \ LU / L iO 



Hr If, \^ , I 2 I 2 I 2 / 2^ 1 , / 2 /2s b 



^'r = , y = 9 1 + , ,2-L X2 (P'rx + Pry + ^'x + ) + T^iPrz + r'z ) + 



The Hamiltonian describing the whole system will be therefore 

Using cylindrical coordinates, that is 

PR = {R^x + R^y)^-, = arctan | ^ j , z/? = i?^. 



/ 2 ^ 



= (r^ + ) " , = arctan ^ = r'^, 

we immediately have the eigenfunctions of H'j^ y given by [39l] 



^J/ijmjjrtjj (r) — _ — R}Jr\mr-\{Pr)s Xn^{Zr) , 
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where R,y\m\{p) are the two-dimensional oscillator radial eigenstates, whose normalized expressions are 

^ { (.Ihy) ' pl"le-«il"l(p^) (43) 

being lL™' the Laguerre polynomials with the quantum numbers v and m taking the values v — 0,1,2,... and 
m = 0, ±1,±2,... respectively. The functions Xr(^) are the eigenstates of the unidimensional harmonic oscillator 
which are given by 

\2"n!7r2 / 

where Hn{z) are the Hermite polynomials and n takes the values n = 0, 1, 2, ... 
The final eigenstates of the Hamiltonian (|38)) will be 

\vRmiinji,Vrmrnr) = \vRmjinji) ® li/rm^nr) (45) 

and the wave function 

(ri,r2) = * (R,r)|J| = {xi,yi,zi\X2,y2,Z2\viimRnii,Vrmrnr) (46) 



where J is the Jacobian of the canonical transformation ((34)) and (1351) . 

The eigenvalues of the harmonic oscillators in one and two dimensions are (n + i) and {2v + |m| + 1), respectively. 
Defining the quantities 



we obtain the eigenvalues of the Hamiltonians H'j^ and H'^ in the form 



KamnnniyR) = ^^^^^ ^ ^ "^^ + 1) + 2^^^^^ + 1"^^' " + 1) + ^^i? + ^ 



(48) 



K.m.njyr-) = y (2^r + |m^| + 771^ + 1) + ^(2i^r + [m^l - + 1) + j^n,. + - j (49) 
Then, the total energy of the system, which is the eigenvalue of the Hamiltonian H, is given by 

E^^mnnn-^^m^nA^^b) = K ^{y r) + (w^ ± \^)'^ E^^^^^^yr) (50) 

We calculate the exact form of the trace of the reduced density matrix associated to a general eigenfunction (|l5)) of 
the two-electron Moshinsky system with magnetic field for the ground and the first excited states in rir, vr and 
Vr- Next we are going to provide and discuss the corresponding amounts of entanglement exhibited by each eigenstate 
(arising from the evaluation of the aforementioned integrals). In what follows, e^,^mrt.nRM,-m^nr denotes value of e 
when evaluated on the state \i>RmjinR,i'r'mrnr) that we also will denote \vR'mRnR,Vrmrnr) Rr- In order to obtain 
physically acceptable solutions in the case of a repulsive interaction between the particles we have to take into account 
the constraint A < 



A. Ground state 1000,000)Rr 

The ground state is symmetric in coordinates, so we must combine it with the only antisymmetric spin function to 
ensure the antisymmetry of the wave function. Let a — — and r = — as before. In this case we have 

8\/TT^\/TT72 VTT^^T^ (2 -t- 2ct2 + r2 - 2%/TT^VTT5^T^) 

£000,000 — 1 , (51) 
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FIG. 3: (Color online) Entanglement of the ground state of the three-dimensional Moshinsky atom with two interacting electrons 
and a magnetic field, a) Attractive interaction, b) Repulsive interaction. All depicted quantities are dimensionless. 



Decoupling the system makes £000,000 = 0, therefore in the decoupled system the ground state is not entangled. 
With maximum coupling t — > 00 (r — > 1) in the attractive (repulsive) case, we find that £000,000 = 1; that is, the 
entanglement measure is maximum. The behaviour of entanglement as a function of the parameters r and a is shown 
in Figure 3. Figures 3a and 3b correspond, respectively, to the attractive and the repulsive cases. More detailed 
information concerning the asymptotic behaviour of entanglement is provided in Figure 4. 

From Fig. ^ it can be observed that in the limit of large magnetic fields, that is cr 00, the entanglement reaches 
a constant value which depends on the relative strength of interaction given by the parameter r, i.e. 



(1 + t2) (2 + r2-2Vr 



lim £000.000 = 1 , , ^ , (52) 



B. First excited states |100, 000)Rr,|000, 100)Rr,|001, 000)Rr,|000, 001)Rr 

We study the excited states in un, Vr and na that have symmetric coordinates wave functions and therefore must 
be combined with the antisymmetric or antiparallel spin function. We also study in this section the excited state in rir 
which is antisymmetric in coordinates and therefore, it can be combined with parallel or antiparallel spin functions. 
Excited eigenstates in u^i and Vr and in and have the same energy respectively when the system is decoupled. 
In this case we obtain 

4\/rT^yTT^VTT^2^^ (8 + 8(7^ + 8r2 + + 8(t2 (2 + r2)) 

£000,100 — £100,000 — 1 ; . a 1 (53) 



£000.001 — £001.000 — t — a n_ , g (,04) 

(1 + VT+^y {VTT^ + Vi + <y^ + T^) 

where a = 1(2) for antiparallel (parallel) spin. 

Taking the decoupled limit system, we obtain the following entanglement values regardless of the magnetic field 

3 1 

value: £100,000 — £000,100 — ^ £001,000 = £000,001 = 2 excited states with antiparallel spin, which are 

entangled, and £gQg = for the only possible state with parallel spin. We have used £" (eP) to indicate the 
entanglement of states with antiparallel (parallel) spin. 

On the other hand, with maximum coupling we find that £?oo,ooo = £000,100 = £001,000 = £000,001 = £000,001 = 1' 
entanglement is maximum. 
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FIG. 4: (Color online) Entanglement of the ground state of the three-dimensional Moshinsky atom with two interacting electrons 
in a uniform magnetic field, a) as a function of r for different values of a with attractive interaction, b) as function of a for 
different values of r with attractive interaction, c) as a function of r for different values of a with repulsive interaction and d) 
as a function of a for different values of r with repulsive interaction. All depicted quantities are dimensionless. 



In the limit for large magnetic fields, i.e. a — >■ oo, we obtain 



lim e 



100,000 



lim £000, 



3\/r 



100 



2(1 + 



2+r2+6Vl+r^ 



2 + T^+6Vl + T^ 

l + yT+72 



(55) 



3 (6 + 3r2 + 2VTTW) J '+^:+/}^ J ^+:l+'ft^ 
lim £001,000 = lim £000,001 = 1 - « , (56) 

which, as we observe, depends on the value of the interaction. The behaviour of the entanglement exhibited by these 
states is shown in Figure 5 (for excited states with Vr = I and vr = 1) and in Figure 6 (for excited states with n,. = 1 
and nji ~ I). 

It transpires from the calculations summarized in Figures 3-6 that the amount of entanglement exhibited by the 
eigenstates of the Moshinsky atom tends to decrease with the strength of the magnetic field. To understand the 
physics behind this trend let us first recall the general way in which entanglement depends on the strength of the 
interaction between the two particles constituting the system. Entanglement tends to increase with the relative 
strength of the interaction. However, it is important to stress that the determining factor here is not the "absolute" 
strength of the interaction, but its strength as compared with the strength of the external confining potential. In 
other words, entanglement increases both if one increases the strength of the interaction keeping constant the external 
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FIG. 5: (Color online) Entanglement of the first excited state in vr and Vr of the three dimensional Moshinsky atom with 
two interacting electrons and magnetic field, a) Attractive interaction, b) Repulsive interaction. All depicted quantities are 
dimensionless. 





FIG. 6: (Color online) Entanglement of the first excited state in ur and rir of the three dimensional Moshinsky atom with two 
interacting electrons and magnetic field, a) Attractive interaction, antiparallel spin, b) Repulsive interaction, antiparallel spin, 
c) Attractive interaction, parallel spins, d) Repulsive interaction, parallel spins. All depicted quantities are dimensionless. 



potential or, alternatively, if the strength of the confining potential is weakened while keeping constant the interaction. 
These general trends have been observed in all the atomic models where entanglement has been studied in detail: the 
Moshinsky model, the Crandall model, and also in Helium and in Helium- like atomic systems [1, [l^. For instance, 
when one considers decreasing values of the nuclear charge Z in Helium-like systems (weakening the Coulombic 
confining potential) the entanglement of the system's ground state increases ^M,- These general patterns admit a 
clear and intuitive physical interpretation. When the external confining potential becomes physically dominant (as 
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compared with the interaction) the behaviour of the system resembles the behaviour of a system of independent, 
non-interacting particles, and entanglement tends to decrease. On the other hand, when the interaction is dominant 
(as compared with the confining potential) the system's behaviour departs from that of a system of non-interacting 
particles and entanglement tends to increase. 

The dependence of entanglement with the magnetic field can now be physically understood. This dependence 
follows the same general patterns explained above. Indeed, the basic fact about the magnetic field in the Moshinsky 
model (which of course is a common external field acting on both particles) determining its effect upon entanglement 
is the following: increasing the strength of the magnetic field tends to increase the confining effect of the combined 
external fields (that is, the harmonic external field and the magnetic field). To illustrate this basic property let us 
briefly consider the behaviour of a single particle (in SD-space) under the combined effects of the external fields 
(harmonic field plus uniform magnetic field) involved in the Moshinsky model that we study here. The probability 
density corresponding to the ground state of the particle is. 



(57) 



A direct way to study the dependence of the confinement of this particle on the strength of the magnetic field is to 
compute the entropy of the spatial probability density and determine its behaviour with the magnetic field (decreasing 
values of the entropy correspond to increasing confinement). The linear entropy, S'*^^-' = 1 — Tr[p^] = 1 — J[p{r)]'^dr, 
and von Neumann entropy, S^""^^ — —Tr[pln{p)] — — J p{r) \n[p{r)]dr , of the probability density ([57]) are given, 

respectively, by ^f^,] [co, 6) = 1 - -fjg+f and sl;^\u;, b) 



3(l-|-ln7r)-lnw-ln (b'^ 



The entropies S, 



(L) 



and sf"^' describing the spatial "spreading" of the probability density associated with ground state wave function 
are plotted against the magnetic field in Figure 7. It can be clearly appreciated that confinement increases with the 
intensity of the magnetic field. The probability densities corresponding to the excited states of this single-particle case, 
as well as the two-particle spatial probability densities corresponding to the eigenfunctions of the Moshinsky system, 
also become more confined when the magnetic field becomes more intense. As already explained, this behaviour is 
consistent with the decrease of entanglement with an increasing magnetic field. 




FIG. 7: Linear and von Neumann entropies as a function of b with cj = 1 



PERTURBATIVE APPROACH 



In this section we consider a perturbative approach to the two previously studied models, regarding the term that 
describe the interaction (between the three or two electrons) as a small perturbation (A^ ~ 0). Let us consider both 
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systems governed by a Hamiltonian of the form 

H ^ Ho + X^H', (58) 

where the unperturbed Hamiltonian Ho takes different forms for each case (see Eqs. ([62|) and ([72l) ). Hq corresponds 
to three independent (non-interacting) particles and two independent particles in a magnetic field, respectively, and 
X^H' describes the interaction between the electrons, being A a small parameter. A perturbative treatment of this 
system involves an expansion of its eigenenergies and eigenstates in terms of powers of . If this approach is valid we 
expect the gross properties of the energy spectrum to be given by the eigenvalues of the unperturbed Hamiltonian. 
It is clear that within this scenario the leading, zeroeth-order contribution to the energy spectrum is independent 
of the detailed structure of the perturbation H' . However, the situation is different when, instead the energy, we 
calculate the entanglement of the system's eigenstates. When the unperturbed energy eigenvalues are degenerate the 
leading (zeroeth-order) contribution to the eigenfunction's entanglement does depend, in general, on the details of 
the perturbation. 

Let us consider an m-fold degenerate energy level of Hq, with an associated set of m orthonormal eigenstates 
j — 1, . . .m. Since Hq describes non-interacting particles, the m eigenstates \tpj) can always be chosen to be 

Slater determinants written in terms of a family of orthonormal single-particle states l^j^'^'"^'') in the case of three 

(1 2) 

particles and in terms of |(/)^ ' ) in the case of two particles. So we have for a three-particle system, 

and for a two-particles systems 

All the members of the subspace Hs spanned by the states are eigenstates of Hq corresponding to the same 
eigenenergy. The different members of this subspace have, in general, different amounts of entanglement. Typically, 
the interaction H' will lift the degeneracy at least partially of the degenerate energy level. If we solve the eigenvalue 
problem corresponding to the (perturbed) Hamiltonian H and take the limit A — )■ 0, the perturbation H' will "choose" 
one particular basis {|?/;fc)A->-o} among the infinite possible basis of Us- The states constituting this special basis will 
in general be entangled. These states are of the form [i^l 



where the m-dimensional vectors — {cki, Ckm) are the eigenvectors of the m x m H matrix with elements given 
by, 

H,j = mH'l^j). (60) 
It is then clear that, in the limit A — the eigenstates of H will in general be entangled. 

("123) (12) 

Let fh be the number of different single-particle states within the family {\(f>j ' ' ), 1, . . . , m} or {j^^ ' ), 1, . . . , m}. 
rh tends to increase with m which, in turn, tends to increase with energy; that is, rh tends to increase as one considers 
higher excited states. This explains (at least in part) why the range of entanglement values available to the eigenstates 
llV'fc) A-fo} tends to increase with energy. Indeed, the amount of entanglement that can be achieved for a given energy 
of a iV-fermion system admits an upper bound given by 



(61) 
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A. Moshinsky model with three electrons 



Let the unperturbed Hamiltonian be, 



19' 19' 1 1 2 2 1 2 2 1 2 2 



Ho = -TTT^TI - oJ^ - o7^ + + n^^2 + o^'^3 (62) 



2dxl 2dxl 2dxl 2 ' 2 ^ 2 

and the perturbation, 



X^H' = \^-[{xi - X2f + (X2 - 0:3)' + {xi - x^Y]. (63) 

Then, we have H = Hq + X^H' . When A = the model consists of three- independent harmonic osciUators with the 
same natural frequency. Let |n±) (n — 0, 1, 2...) the eigenstate of each of these oscillators. For the first excited state 
which is four-fold degenerate, let {|0, ±), |1, ±), |2, ±)} be the single-particle orthonormal basis. Then, for A = 0, 
we can choose the four eigenstates with zero entanglement, all of them with the same energy as \011) li^ji^R:^ and 

|llO)_Rifl2-R3 

iV'i) - |0+,0-,2-H| 
|V'2> - |0+,0-,2-| 
|V3> = |0-H,1+,1-| 
1^4) = |0-,1+,1-| 

(64) 

where we have introduced the notation 

\hj,k\ = {\hj,k) - \i,k,j) + \j,k,i) - \j,i,k) + \k,i,j) - \k,j,i)) 

and i = J = fc = 0±, 1±, 2±. 

For the first excited energy level of Ho [E = we have 



H cx 



/ 4 
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V2 
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7 
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1 


7 

2 y 



(65) 



and the corresponding eigenvectors can be written as 



K) = y|(-^iV'2) + iv4)) 

W2) = /f(-^l^l) + l^3)) 

1^^) = -1(721^2) + 1^4)) 

W,) = -L(y2|^i) + |^3)). (66) 

In the decoupled limit the eigenstates r^r^Rs ^"^^ \^^^) R1R2R3 tend to the states (or any combination of 
them) which have e = ^ • This value coincides with the amount of entanglement for the first excited state obtained 
from the exact calculation in the limit r — )• (equivalently A — )■ 0). 

The states |V'i) and IV'2) (IV'3) and IV'4)) share the same energy eigenvalue. A linear combination of eigenstates 
sharing the same eigenenergy is also a valid energy eigenstate. Then, let us consider for instance IV'34) = cosO \^'^) + 
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sin Olip'^}, {0 < 6 < 2tt). As already discussed at the end of Subsection III.B, the amount of entanglement of these 
linear combinations does not depend on 0. 

We have seen that in the case of some excited states of the Moshinsky model an arbitrarily weak interaction between 
the particles leads to a finite amount of entanglement. This naturally suggests the following issues: to what extent is 
this weak-interaction entanglement robust? What happens with this entanglement if some other small perturbation 
acts upon the system? The detailed entanglement features of the eigenstates corresponding to this scenario will 
evidently depend on the precise form of the new perturbation. Therefore, these entanglement properties can only 
be studied in a case-by-case way. However, it is possible to gain some valuable insights on the robustness of the 
weak-interaction entanglement by recourse to a statistical approach. We can consider the typical features of the 
weak-interaction entanglement corresponding to a random perturbation. 

Let us consider again the four eigenstates (1641) of the unperturbed (with no interaction) system. Any weak pertur- 
bation acting on top of the already considered weak interaction will lead (in the lowest order of perturbation theory 
for a degenerate eigenenergy) to a new set of four perturbed energy eigenstates that will be linear combinations of 
the four unperturbed states (l64|) . That is, the new perturbed states are orthonormal states belonging to the four- 
dimensional linear space spanned by the states (j64| . We can consider the statistical distribution of entanglement 
values corresponding to random states in this subspace uniformly distributed according to the Haar measure (see 
[111,113 ^'^d references therein). To this end we generate three-electron states randomly distributed according to the 
Haar measure of the form, 



|V'') = 5^q|V'.), (67) 

with \tpi), i = 1, . . . , 4 as given in Eq. ([M)) . A state of the form \tp') can be thought as an eigenvector corresponding 
to an arbitrary perturbation. The amount of entanglement of the state \ip') is. 



H\ci\ 



|C2| 



2(|C3| 



|C4| 



(68) 



Optimizing Eq. (|68p we obtain the maximum possible value of entanglement EmdV'')) — ^ associated to the state 
with coefficients satisfying: |cip -I- |c2p = |c3p -I- |c4p = ^■ 

In Table U we show the percentual number of three-electron pure states belonging to the linear subspace spanned by 
the states (j64p that have entanglement values in different ranges. To compute these percentual values we generated 
10^ random states distributed according to the Haar measure. The average and maximum entanglement values 
corresponding to perturbed states spanned by (|64l) are also given in Table HI These results constitute suggestive 
evidence for the robustness of the entanglement exhibited by excited eigenstates of the Moshinsky atom in the weak- 
interaction limit. Indeed, the statistical study summarized in Table |T] suggests that any new perturbation is likely to 
produce a small decrease in the entanglement of the excited state considered here (e = 8/27 w 0.2963) but, in the 
typical case, it will still result in an appreciable amount of entanglement in the weak-interaction limit. 



% of states in 
different 
entanglement ranges 


< e < 1/9 


4.75% 


1/9 < e < 2/9 


18.25% 


2/9 < e < 1/3 


77% 


average entanglement 


(e) = 0.26667 




maximum entanglement 


Em = 1/3 





TABLE I: Entanglement distribution for perturbed excited states of a three-electron Moshinsky system. These states correspond 
to a four-fold degenerate unperturbed energy level. 



Let us now study the entanglement properties for the second excited state of the Moshinsky atom. In this case, we 
have ten- fold degenerate eigenstates, (all of them with E = |a;). The single-particle orthonormal basis is given by 
{|0,±), |1,±), |2,±), |3,±)}. The matrix of the harmonic perturbation can be expressed as follow: 
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(69) 



Following a similar procedure and after a laborious algebra we compute the entanglement amount of the eigenvectors 

of 



e{W,))^s{W,))=Q 
20 



(70) 



where \ip'j) [j = 1, 10) are the eigenvectors of the H matrix. The states \ip'j) with j = 1,...,6 share the same 
eigenvalue. The same occurs for the state pairs lipj) and \ip'g}, and Itpg) and iV'io)- we mention before, the 
interaction lift only partially the degeneracy. The degeneracy due to the spin degree of freedom {Sz = ±^) is present 
in all the states. The obtained values agree with some of those calculated in section 3. Different combinations 
of the states sharing eigenvalues result in the non-coincident entanglement amount eo2i — For instance, let 

iV'^e) = pVs + Vl-p'V-e, (0 < P < 1); then 



(71) 



and £021 is re-obtained for p ^ 0.992. 



B. Moshinsky model with two electrons in a uniform magnetic field 

We consider also a perturbative approach for a three-dimensional Moshinsky atom with two electrons in a magnetic 
field. Let the unperturbed Hamiltonian be. 



y (a;? +yl + xl+ yl) + h{Liz + La^) 



and the perturbation. 



(72) 



(73) 



The eigenenergies of are given by Eq. (jSOp . taking vji — vi, Vr — 1^2, ttt-r = mi, nir = m2, — ni, rir = n2 and 
setting A = 0. Then, for the excited states of Hq with energy given by 



£:^,„=a;(2y+^ + l^ , y = j^l + + A, 



19 



resulting of setting one of the quantum numbers ui, 1^2, \m2\ equal to one and the rest equal to zero, we obtain 
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where ci = ^ + ^^^^^^ and C2 is obtained numerically and its exact numerical value is not relevant for the next 
calculations. 

H has six degenerate eigenvectors and two non-degenerate ones that take the following entanglement values: 
{0, ^,|}. The entanglement value obtained from the exact computations in the limit A — >■ of the states with 
the same energy, |100,000)flr and |000, 100)/i'r5 coincide with one of the above values (e = |). 

We consider also states setting rii = 1 or 712 = 1 and the rest equal to zero, with energy given by 

and with y as before. For these excited states we obtained 



H<x 



(di+d2 



V 









di 


d2 


d2 


di 
















where di = — 

J- (.1 



and d2 = This matrix has three eigenvectors corresponding to the same eigenvalues and 



with entanglement {0, ^} and one non-degenerate eingevector with entanglement i. Again the obtained results are in 
perfect accordance with the entanglement obtained for states 1001,000)^^ and |000, 001)i{r in the decoupled regime. 



VI. CONCLUSIONS 



We explored the entanglement properties of two versions of the Moshinsky model: one comprising three electrons 
and another one consisting of two electrons in a uniform external magnetic field. The amount of entanglement of the 
eigenstates of the three-electron system considered here depends only on the dimensionless parameter r describing 
the relative strength of the interaction between the particles (as compared with the strength of the external confining 
potential). We obtained closed analytical expressions for the amount of entanglement of the ground, first and second 
excited states. As a general trend we found that the entanglement exhibited for these states tends to increase 
both with the state's energy and with the strength of the interaction between the particles (that is, with r). Non- 
vanishing entanglement is obtained in the limit of vanishing interaction in the case of excited states. This (apparent) 
discontinuous behaviour of the entanglement is related to the degeneracy of the energy levels of the "unperturbed" 
Hamiltonian describing non-interacting particles. The non-vanishing entanglement in the limit of zero interaction is 
determined by the particular basis of Hq "chosen" by the interaction. We also found that in the case of an attractive 
interaction the eigenstates' entanglement approaches its maximum possible value in the limit of an infinitely large 
interaction. On the other hand, in the case of a repulsive interaction the maximum possible entanglement is obtained 
when the interaction strength approaches a finite, critical limit value corresponding to Tc = ^ . The system does not 
admit bound eigenstates when the strength of the (repulsive) interaction is equal or larger than the one corresponding 

to Tc. 

As far as the entanglement's dependence on the interaction strength and the energy are concerned, the behavior 
of the Moshinsky model with two electrons in a uniform magnetic field is similar to the one observed in the three- 
electron model. With regards to the external magnetic field, we found that the eigenstates' entanglement decreases 
when considering increasing magnetic fields. In the limit of very strong magnetic fields the entanglement approaches 
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a finite asymptotic value that depends on the interaction strength. The essential aspect of the magnetic field in the 
Moshinsky model that determines its effect upon the amount of entanglement exhibited by the system's eigenstates 
is the following: increasing the intensity of the magnetic field tends to increase the confining effect of the combined 
external fields (that is, the harmonic external field and the magnetic field). For a given strength of the interaction 
between the particles, this increasing confinement leads (according to a general pattern that has been observed in all 
atomic models where entanglement was studied in dcitail) to a decrease in the eigenstates' entanglement. As happens 
in the case of the three-electron model, a perturbative treatment highlights the essential role played by the degeneracy 
of the energy levels of the interactionless system in determining how the eigenstates' entanglement depends on the 
interaction strength and on the energy. 
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